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Abstract 

Off-shellness and inhomogeneities in the non-equilibrium dynamics of the A(/3^ model are studied. The 
non equilibrium field theory is developed in the Closed Time Path formalism, which is in turn reformulated 
as a kinetic field theory, in terms of a set of kinetic equations for the 2-point Green function, derived from a 
generalized effective action functional. We take into account all initial correlations up to the 4-point functions. 
It is shown that the model displays an SO(l,l) symmetry broken by interactions and initial conditions. The 
divergence of the corresponding Noether current is identified as a particle creation/annihilation density. The 
broken Ward-Takahashi relations for the SO(l,l) symmetry are derived. They constitute a set of generalized 
kinetic equations for the general n-point functions. Energy-momentum conservation is demonstrated to hold 
from the fully interacting transport equations, and wc discuss the effect of a gradient expansion and its non- 
locality. As an application, we study an explicit realization of inhomogeneities and non-equilibrium conditions 
in the free field model, for which we give the general solution and analyze it in a particular case. The 
identification of the Casimir effect in our solution illustrates the intimate connection between inhomogeneities 
and off-shellness. 
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1 Introduction 



Phase transitions play an important role in the study of physical systems under extreme conditions. In 
particle physics, various kinds of phase transitions are of interest and the use of equilibrium thermody- 
namics may be inadequate. For example, this is so in the current descriptions of the initial stages of heavy 
ion collisions, in which, heuristically, the onset of deconfinement is believed to give rise to the formation 
of a quark gluon plasma [Q. Two further examples that one may speculate upon, are the electrowcak 
phase transition in the early universe [|| , and the formation of disaligned chiral condensates (dec) in high 
energy hadronic collisions [Q , such as tt p scattering. 

The closed timepath (CTP) formalism, originating from the work of Schwinger 0, Keldysh ||], and 
Vernon and Feynman is suitable for analyzing non-equilibrium systems. It takes one beyond the more 
conventional approaches used to describe small perturbations about thermal equilibrium, such as thermo 



field dynamics Q, the Niemi-Semenoff methods |10| and the real time finite temperature formalism of 
Dolan and Jackiw The CTP formalism is not restricted to quasi homogeneous systems close to 

equilibrium, although in practice, to date, it has been difficult to implement non-equilibrium theory in 
inhomogeneous systems. The full power of the CTP formalism is unfolded within the functional integral 
representation ||l2|, |l^ . Recently |l^ an effective action principle, that leads to a hierarchy of Schwinger- 
Dyson equations for the Green functions, has been formulated. These equations can be reformulated in 
terms of a transport theory, or kinetic field theory, for the 2-point functions. Such a field theoretic 
formalism can be shown to have a Boltzmann equation as its quasi-classical limit While the CTP 
non equilibrium formalism is theoretically aesthetic and complete, it is difficult to implement practically, 
and has only recently become a subject of detailed study. Transport theories in general constitute a topic 
whose quantum field theoretic content is difficult to understand on the basis of quasi-classical models. 
We briefly summarize some applications of transport theory that have been developed and are currently 
of interest. We do so without attempting to be complete. 

Non-relativistic kinetic equations have become a common tool in analyzing heavy ion collisions nu- 
merically p^ . Semi-classical, but relativistic, approaches based on quantum hadron dynamics have been 
put forward for the same purpose |^ . A lot of work has been done on the formulation of a quark-gluon 
transport theory based on the operator formalism of quantum flcld theory ]l8| . In order to study the 
restoration of chiral symmetry in relativistic heavy ion collisions, a transport theory for the Nambu- 
Jona-Lasinio model, that contains quarks as its fundamental degrees of freedom, has been derived in the 
framework of the CTP formalism [Q. Due to the difficulties inherent in this particular problem the 
equations of motion that are obtained have however only been solved in the mean field approximation, in 
order to study the influence of the phase transition on the expansion of a QGP and the fluctuations 
in a quark-antiquark plasma |2l[. The CTP formalism has also been recently applied to the 0(4) 
linear sigma model to investigate dec formation. In all of the studies that have been mentioned, the sys- 
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terns have been required to be quasi uniform and/or the particles to be on mass-shelL These restrictions 
have simply been necessitated by the dificulty of the problem at hand, since they give rise to considerable 
simplifications in the formalism when applied. On the other hand, the lifting of these requirements may 
in some cases be important. The mass-shell restriction may suppress important quantum effects. This is 
illustrated in the calculation of Eisenberg jl^ showing that the spontaneous creation of particles from a 
constant background field cannot be described in this way. 

In order to obtain a better understanding of the importance of off-shell behaviour and spatial inho- 
mogeneities, it is useful to examine a simple model. To this end, we shall focus on the well-known scalar 
Xif^ model in this paper and study the CTP formalism and associated transport theory for this model. 
In doing so we shall endeavour to remove the on-shell restriction, and also investigate a scenario which 
is far from equilibrium. Of technical importance is the proof that energy and momentum are conserved 
in the presence of interactions, which we also demonstrate. 

Since our analysis is performed in the setting of kinetic theory, we discuss this briefly here. If inter- 
actions are short ranged, as is the case for the Xip"^ point-like interaction, one expects naively that the 
quantum field may, above some length-scale be conceived as a distribution of free particles, each of 
which temporarily goes off-shell while undergoing an interaction, during a time . For this picture to 
have some validity, one has 

Tc«^/. (1.1) 

In that case, the time dependent particle density n{x, t) satisfies 

n(f,t)<^-3 (1.2) 

in describing the system. Due to particle creation and annihilation n[x, t) will, in general, have a time 
dependent normalization. 

This separation of the system into macroscopic length scales above € / , and microscopic length scales 
below is at the same time a statement about the effective interaction, and hence the effective coupling 
of the particles at different energy scales. It is therefore important to consider the behavior of the 
running coupling as we move from one length scale to another. For example, in QCD at high energies, 
the quark and giuon degrees of freedom are relevant, whereas at low energies these remain confined 
and their hadronic bound states can be treated as fundamental. In this case, the separation in length 
scales is accompanied by a separation of physical degrees of freedom. Such a phenomenon is difficult to 
implement in a two-scale transport theory. A two-scale approach is also not consistent with the running of 
the effective coupling in the Xip"^ model. Here we are dealing with a trivial theory. If we take the quantum 
corrections of all scales into account, the effective coupling in the theory vanishes, implying that small 
scale dynamics fundamentally modify the behavior at larger scales. Imposing a momentum space cut 
off in the theory invalidates triviality and the theory becomes interacting. A kinetic formalism can be 
derived from the closed time path formalism without resorting to a two-scale interpretation. However 
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the resulting equations of motion that one obtains are kinetic equations for the 2 point Green functions, 
rather than for a quasi-classical phase-space distribution function. 

A kinetic theory for the Xip'^ model was originally derived via the operator formalism in the mid- 
seventies [Q . Wc study the kinetic field theory for the aforementioned model by reformulating the CTP 
equations of motion as a transport theory, and including up to 4-point initial correlations. A similar 
approach can be found in where however only initial two-point correlations have been taken into 
account. In section two, we give a brief introduction to the CTP formalism in this context, derive the 
transport equations, and clarify their physical content. Our view parallels the physical interpretation 
put forward in the early work by Carruthers and Zachariasen Initial conditions for these equations 
are shown to break a symmetry underlying any CTP field theory in the third section. The Ward- 
Takahashi identities corresponding to this broken symmetry express direct relationships between the 
Green functions and the initial conditions. They also represent a set of kinetic equations generalized to 
the n-point Green functions. For the A(p^ model, in which the symmetry group is S0(1,1), these identities 
are given. In particular, the Noether current corresponding to this symmetry is interpreted physically. In 
addition, it is shown that the dynamical evolution of the 2-point functions guarantees energy-momentum 
conservation. In the fourth section we detail this explicitly by demonstrating that the divergence of the 
energy- momentum tensor vanishes simply by using the transport equations. This allows us to settle 
some of the questions concerning energy-momentum conservation and the use of gradient-expansions in 
both Vlasov equations, and equations that include a collision-term. By studying the off-shell behavior in 
the free field equations and recovering the Casimir-effect from our solutions, we illustrate the intimate 
connection between inhomogeneities, off-shellness and initial conditions. With this program in mind, we 
now give a brief exposition of the CTP formalism and the kinetic field theory that is based upon it. 

2 Kinetic Field Theory 

Before we turn to the Xip^ model in particular, we review the general structure of the CTP formalism 
and the generalized effective action. Following this, we derive the transport equations for the 2-point 
functions and analyse the physical meaning of the various terms in the equations and the variables we 
use. 

2.1 The Closed Time Path FormaHsm 

Here we will briefly recapitulate the main ideas of the CTP formalism that are relevant for our calculation. 
A full treatment of the CTP formalism can be found elsewhere Let (p denote a collection of field 

operators. For simplicity, we take them to be bosonic. The system is specified by a Lagrangian density 
C{(p), and a density operator p, that enforces the initial conditions. 
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Let I 01 > be an eigenstate of the ip operator satisfying 

^> |0i=<^i(a;) \(j3i> (2.1) 

with the c-number function (j^iix) being one of its eigenvalues. Let | 0„ > be a complete set. Then the 
functional 

p[Ux),^M] = < I p(t = 0) I > (2.2) 

contains all the physical information stored in the density operator at t=0. Square brackets denote p 
being a functional rather than a funtion. If the hamiltonian TL of the system is not time dependent then, 
in the Heisenberg picture, neither will the operator p be. If we now introduce an additional interaction 
with a source j{x, t) 

H^n + J d^xj{x,t)^{x,t) , (2.3) 
then the time dependence of p is given by 

p{t) = [T{Gxp(-z J'dt' J d^x'j{x',t')ip{x',t') )} ] 
X p{t = 0) X 

[r{exp(+* /* dt" I d^x"]{x'\ t")ip{x", t") )} ] (2.4) 







where T is the time ordering operator, and T is the anti time ordering operator. It is a useful bookkeeping 
device to introduce indices 1 and 2 for the current and field operators taken on the time ordered or anti 



time ordered paths respectively, that are given in Eq.(2.4). Explicitly on the time ordering path, one 
has (p^(x,t) and j^, and ifP'{x,t) and on the anti time ordering path. In a contravariant notation, the 
general indexed n-point Green functions are defined by 

e"-'^"(a;i,...,a;„) = {Tpi^'^' {x^)...^'^- {xn)}) (2.5) 

where Tp orders the field-operators according to the rule 

%{\{ v"- M } = ni[ {xn)}T{ n i^k)}. (2.6) 

n an— 2 0.^^ — 1 

The action = J d'^x C enters in the functional integral representation of the trace of the density 



matrix, that is obtained by writing the time-evolution operators in Eq.(2.4) as path- integrals, and using 



the functional Eq.(2.2) 



2\j\p] = Tr{pit = ^)} 

= J D^' D<j>^ exp + o 0i - - f o ^ (2.7) 

using the abbreviation j o 4> — J d'^xj(x)(l){x). 



4 



It is convenient to replace the fields (j>^ and 0^ by the doublet (p°- — {(p^, (fp'), and introduce the metric 

Cab = c''^ = I ^ M (2.8) 
on this internal space. For example, the source term can then be written as 

jl ° 01 - i2 °(t>2= Cabf OLp'' = o ip" (2.9) 

by contracting the fields (j) with the metric c°^. We will refer to these indices a,b as CTP indices. We 
define a CTP action Sctp = S[(j)i] ~ >S'[02] and write for the partition function Z 

Z[ja, p]^ J Dip' exp liScTpif'] + Ja o }P[A ■ (2.10) 

The general Green functions Q'^'"^- - are computed now by taking derivatives of Z with respect to the 
source ja (x) . 



By expanding the density functional in Eq.(2.2), in powers of the field (p" 

oo ^ 

p[ip^] = exp -,Ka,...a^ O O ... O } , (2.11) 



n! 

ri=0 



one obtains the expression for the generating functional of the general n-point Green functions CJ'^i - '^" 
to be, 

Z[ja,Kab,-]= J D^''exp l{ScTp[V^'']+jaOv'' + f2^Ka,...a„0^''«O...O^aJ. (2.12) 



From Eg. ( 2. 12 ), one sees that the non-local sources Kabc... act as sources for the full Green functions g'^'"^ --. 
In this way they prescribe the initial values of these Green functions. In a diagrammatic expansion of 
Eq.( ^.12" ), the non- vanishing K's are part of the tree level contributions to the non-equilibrium vertex 
functions, although they contain contributions from all loop orders. 

The generating functional W of the connected Green functions G"''^ ' (a;, y, z, ...) is the logarithm of 

Z 

W[ja,Kab,...] = -t\0g{Z[ja,Kab,...]} • (2.13) 

and the connected Green functions can be derived from this in the standard fashion. An alternative 
method of obtaining the connected Green fuctions is as follows. One derives the full Green functions 
Qabc... via differentiation of W with respect to the non-local sources 

g-^' -^--{xu..:X„)=n\— U,if6c,...] , (2.14) 

and obtains the connected Green functions by subtracting out the disconnected contributions. This 
procedure is most simply convenient in deriving the hierarchy of coupled equations of motion for the 
(connected) Green functions in a direct fashion. In order to obtain these equations for the G'^''^", and 
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Qab particular, we construct an apropriate action functional by means of a Legendre transform with 
respect to the non-local sources. Let us denote the mean- field G°'{x) as (^^ . Defining Fas the Legendre 
transform of W, 

n^'^ {x),G''\x,y),...]^W[3a.Kab...]-jaO (p- - l^f^, o {G"^" -|- (p- ^"1-... (2.15) 

the mean-field (p"" and the connected 2-point and n-point functions satisfy the equations 

(5r 



= -ja-KabO - (2.16) 



n\ - ■■ SG'^''' 

n=2 



n—m 



(2.18) 



The general Green functions Q"^'"^- - can be expressed in terms of the connected Green functions G°^^'^^ ^ 
yielding the formal equations of motion. They are formal in the sense that F has not been computed 
explicitly. For our purpose this will not be neccesary. The physical implications of such a Legendre 
transform were treated, for example, in JTzt and [ p6| . At this point we leave the formal development of 
the CTP-formalism and turn to its application to the Xip'^ model. 
We consider the Lagrangian density 

C = \d,^d^^~'^^^^^V^ . (2.19) 



The corresponding CTP action is then 



ScTP^'^] - j d^x ^Cabid^^^'d^^^" ~ mW"} - ^ Wf^WV" (2.20) 
where we have defined 

{1 when a^h = c = d= \ 

) ■ (2.21) 
— 1 when a = b = c = d = 2 I 

In the next subsection we focus on the equations of motion for the 2-point connected Green functions, 

since they lie at the heart of the kinetic formulation of the non-equilibrium theory. 

2.2 The Kinetic Formulation 

In this subsection we derive the equations of motion specifically for the 2-point functions G°^(a;,y) and 
write them as kinetic equations. In particular, a hierarchy of coupled equations for the G"''^' ■ is obtained 
in which the two-point function G°^ is embedded. For practical calculations, this hierarchy of equations 
for the G°^'^ " must be truncated. Calzetta and Hu §, Kandrup and Hu ||2^ and, in a general context, 
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Hu p5[ , have argued that this introduces dissipative phenomena in the modeL On the other hand, if 
the field is initially free and the density operator is diagonal in momentum space, an explicit calculation 
of the non local sources K is feasable, and shows that the series truncates after the two point source. 
Consequently if one is willing to accept limitations on the range of initial conditions, the series of the K's 
can be truncated. It is in this spirit that we will truncate the expansion of after Kabcd{u, v, w, x). 

2.2.1 The Effective Action for the Ai^" Model 
Truncating F after the fourth term, leads to the explicit form 

r[(^'^ ,G"'^G""'^G'•^^'*] = mJa,Kab,Kabc,Kabcd] (2.22) 

'JaO ^'^ -^KaboG^' ^j^KabcOG^'-- ^KabcdoG"'"' 

implying that , G''^{x,y), G''^''[x,y,z) and G''^'"^{w, X, y, z) can be chosen as independent variables. 
A more convenient choice are the 2-point connected Green functions G°''^{x, y) 

G''\x, y) = 2G''\x, y) - (x) (y), (2.23) 

and the 3- and 4- irreducible vertices of the theory, a^(,j,(a;, y, z) and ot^ijcdi'^^ ^) respectively, defined 

by 

ial^^oG'"''G'''''G'"=' ^ilG"''''''' -G"'''' ifi'' -G^'"' ~ G"' >f''' - (p"' (p^' (p"' (2.24) 

and, 

Jaabcd ° G""' G''''' G""' G"'''' + G""' G^^' o ial^„^ o G™"' o la^,^^ o G""' G'^'^' + permutations 

= 4\g"^'b'c'd' _ disconnected terms , (2.25) 

and which are represented diagrammatically in fig.l . As shown by Weldon p8[ |, in the non-equilibrium 
formulation these vertices are related to the equilibration rate of the distribution of quasi-particles over 
the available energy-momentum space, and not to the decay-rate for a single particle. 

2.2.2 The Equations of Motion 

It is possible to extract all information needed to construct the equations of motion for G°^ from the 
dynamics of the mean field ip'^ alone. For the sake of convenience we will always assume that the sources 
Kabc and Kabcd include the tree-level vertices. We proceed as follows. The equation of motion for is 
found by solving the identity 

f S ^ 1 

= D^'' ^exp i{So CTP[/] +jaO^'' + Y, ^Kao...a„ ° ° - ° v""}, (2-26) 
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where the free CTP action 5*0 ctp, given by 

So CTP = J d^x^^id^^-d'^v" + mW'}, (2.27) 
appears due to our choice of absorbing the interaction vertices in the K's. Working this out yields 

Gro) a, - = 3a + E ;jy^aa....a„ O (2.28) 



where we have denoted as G^^q^ the inverse of the free propagator 



From Eq.|2.15|, it follows formally that 



93 2 (5 i^" 3! 



Comparing Eq.| 2.28| and Eq.[ 2.30 |, and using 

§gaa2...a„ 



we can identify the l.h.s. of Eq.|2.3C] to be 

ST 



(2.29) 



(2.30) 



(2.31) 



(2.32) 



^(0) ab "^^^ ^ functional of the full propagator. The explicit connection is given by the usual Dyson 
equation. As in the zero-temperature formalism, Eq.[2.32| can be formally integrated with respect to 
ip"" yielding 

r[ (p- , G''^ o^lbcd] = ro[ , G"^"] + ^.„i[G"^ o^tbcd] ■ (2.33) 

Tint can now be constructed from vacuum diagrams alone, in terms of G"^, a^i,c ^^'^ '^abcd ' which 
simplifies the calculation considerably. To see this, consider, for example, a diagram D pertaining to Tint 
with propagators, 3-point vertices and n4 4-point vertices. It has no external lines, since they 
would end in an cp"" , so 

^ (2.34) 



no = 2'*^'^"^ '^^'^^ ■ 

Let Vint{D) be the Feynman-amplitude of the diagram D contributing to the functional Tint- Cutting 
through one propagator line in D represents the action of a functional derivative, with respect to G°^, on 
Tint{D). The operator G"'' o j^tk^ thus counts the number of propagators in every diagram. Using the 



appropriate counting operators for the vertices, and using the fact that Eq.|2.34| holds for all diagrams 
contributing to Tint, one finds 



Qab 



^ P 33 

o -rrr-^l int = -a 



r,^abc r 3 Tint + 2a j,^ < 

2 Sa'lt,^ 



5at 



(2.35) 



bed 
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The formal equation of motion for G , after truncating the series of the Kabc, is 

_ ^ 3 4 1 1 1 7^ ^Gabc 1 ^Gabcd ,„ oa\ 

-g^l V , a^fc^, a^bcd] = "2 ~ gj^abc o - ^i^abcd o -^^j;^ • (2.36) 

The l.h.s. of this equation can be expressed as 

^[ ,G^^aL,«LJ - ^[ ,G=^«L,«'aU + ^[ (2.37) 
The first term here can be exphcitly calculated from the functional integral for Z 

ro[ al,,, al,J = So ctp[ ^''\~\ MG''''} + U}-^^ o (2.38) 
The last term can be written out using the formal equations of motion of the vertices a^i,^ 

= _1G"'^G''^G^/ o {Kdef - Kdefg 0' } ~ ^G'^'^G^^ O K^efg o Gf^^G^' o a^, ° (2.39) 



ahc 
4 



^ = -G^^G^^G^f G'''' o K.fgh, (2.40) 



obtained from Eq. |2.15|. After some tedious algebra, we find the equation for G°^ to be 

^G-J = ^G^o5 - habcd o { + G^-^} 

-#WdO a^/.^.oG^^GS^G''^ 

-^^hbefgo ^SaL<j°G-G'^^ (2.41) 

- A/i,,,, o al f^a%^ o G^^G^'G^'^G''^ 
-|/i,,,e° o G/^G'"'G^"^G''^ 

The solution of this equation requires knowledge of both vertices a^^^ and a^jj.^;. It describes the exact 2- 
point Green function for the non-equilibrium system, with initial conditions constrained by the truncation 
of the K's after Kabcd- At this point, it makes no explicit reference to the dynamical role of inhomogeneities 



and off-shellness. The kinetic representation of Eq.[ 2.41 however allows one to address these questions 
directly. Since we aim at understanding the simplest effects of inhomogeneities in the field theoretic 
setting, there is no apparent need for allowing spontaneous symmetry breaking. This allows one to make 



an appreciable technical simplification of Eq.[ 2.41 ]. We thus set the mean- field and the 3- vertex to zero, 
so that the last three terms vanish. 



We progress toward the kinetic representation of the non-equilibrium dynamics by multiplying Eq.(2.41 
once from the right by G^'^ and contracting over the index h yielding 

l5,^ = zG(o5 ab ° - \habef O G'^ G'^ - '-^held O aafgh o G^' G^'' G^"' G'' . (2.42) 

Here the differential operator Gj^j acts on G'^''. Multiplication of Eq.( 2.41 ) from the left by G^° and 
contracting over the index a yields 

= tG-- o G^„i - ^habef o G^-G'f - '-^G'^^hbeid o aafgh o Gf'GS'^G^-. (2.43) 
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In doing so we have formally inverted the operator G^i^{x,y). In order to do so properly it is neccesary 
to specify its kernel, i.e. the set of functions f{x) for which 

Gj{x,y)of{x,z) = 0. (2.44) 

This coincides with specifying boundary conditions on Eq.( p.42| ) and Eq.(2.43). Allowing inhomogeneities 
in the boundary conditions entails the need for both Eq.(2.42) and Eq.( 2.43| ) instead of simply Eq.(2.41). 

For a homogeneous system , G°^^{x, y) depends only on the relative coordinate x = x — y. It is there- 
fore natural to consider G"'' as a function of x and the centre of mass coordinate X — ^{x + y}. Its 
4- Fourier transform G°-'^{P, p) 

G''\P, p) = j d^Xd^ X e'^^+' ^ PG^\X, x) 

is related to the Wigner transform G°-^{X, p) by 

G'^\X, p)^l ^^e-^P^G^\P, p) . 

Then Eq.(^.42|) and Eq.(^.43[) in terms of G°''(P, p) become 

d^P' d^ p' 1 



<P+f)'G"''(P, p) 



A 
2 

12 7 (27r)^ 



{27r)4 (2^)4 
d^P' 



G'^'^iP', p')eM^P'-dp}G'^\P-P', P) 
T.\{P-P', p)A-G^\P', p), 



and 



,i^p.PfG^\P, p) 



,')exp{-V.a,}G-(P-P', p) 



A 

2 _ 

i\ r d'^p' 

127 (2^ 



S]",(P-P', p)A+G"''(P', p) 



(2.45) 



(2.46) 



(2.47) 



(2.48) 



where the 2-loop self-energy, T,ab, is given by 



hbcde o al,,,,,, o G- G^^ G 



(2.49) 



and we have introduced the gradient operators 

1 



A±=exp{±-{(P-P').a^-P'.9i} 



(2.50) 



2^^- - p --pi 

In this expression, the arrows indicate the direction in which the differentiation acts. By subtracting 



Eq.(2.43) from Eq.(2.42) we obtain the kinetic equation 



I p.PG^'iP, p) 



A 
2 



WfWf^'^'^^'^'' P>i^M^^'.5p}G'^^(P-P', p) 
d^P' 



12 J (27r)' 



{^''^A-G'''' -^J'A+G"''} 



(2.51) 
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and, by averaging the equations, Eq.(2.43) and Eq.(2.42), the constraint equation 

A f d^P' d^p' 1 



2 J (27r)4 (27r)4 



12 J (27r)4 
ic''\2nf5^{P) , 



G""'(P', p') cosh{- P' .d p}G''\P - P' , p) 



(2.52) 



is obtained. The so-caUed gradient expansion for the Green functions can now easily be obtained by 



a truncation of the series expansions of the operators in Eq.|2.50| or alternatively the sinh and cosh 



occuring in the Eqs.|2.51| and |2.52|. It is a well-defined approximation whenever | P^d^p \ is small. Since 



P' is integrated over, this suggests that the P'-integration is regulated at an ultraviolet limit A and the 
gradient is restricted similarly, i.e. \ d p \< A^^. Such a cutoff neccesarily introduces a length scale. 
Furthermore it also alters the nature of the Xlp'^ model substantially. Both topics have been discussed 
in the introduction. If the operators are to act on functions containing resonance-like singularities in p, 
then the approximation is no longer valid close to the pole for any value of P. 

We close this section with a discussion of the physical content of the various terms in the kinetic and 
constraint equations. 



2.2.3 Physical Interpretation 

Under homogeneous boundary conditions, such as global thermal equilibrium (T.E.), we have 

G'^\P, p) (X 5\P) , (2.53) 



which, when substituted into Eq.(2.51), makes the kinetic equation void. The constraint equation reduces 
to the Dyson equation. The variable p can then be identified with the momentum of the propagating 
boson, and, for example, G^^( p) can be written in terms of a Bose-Einstein distribution function /( p) 
as 



G^\p)^{e{p^) + f{pmp''-m^) . 



(2.54) 



Similarly, one can identify p.P G {P, p) as the flow-term of the kinetic equation. The terms containing 



the sinh function in Eq.(2.51) describe the effect of the interactions with a background of spontaneously 
created (off-shell) particles. Due to the analogy with the Vlasov interaction in plasma physics, we refer 



to it as the Vlasov term. The Vlasov term of Eq.(2.51) describes a strictly local interaction of a single 
particle with the density of other particles. There is no long range component to this. This situation 
may change however when a mean field (p"" is present that contains also non-local information. 

The term in the r.h.s. of the kinetic equation is the collision term. On-shell, it describes the effect of 
two-particle interactions on the one-particle distribution function /. The relation between the two- loop 
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self-energy on the one hand and the decay of the bosons is found from the cutting equations implementing 
unitarity 63] . The Boltzmann equation for a Bose-gas is obtained from the kinetic equation by inserting 



the form of the G"'' from Eq.|2.54| into Eq.(2.51), generalizing /( p) to f{P, p) and neglecting all gradients 
in p. This equation then describes on-shell bosons. The previous remarks on the validity of the gradient 
expansion apply here. On-shellness is not compatible with neglecting gradients. The problem can be 
softened by asssuming that the boson self-energy has an imaginary part, so that ||2^ the distribution / 
is not in equilibrium. In that case, the S{ p^ — "w?) would be replaced by a Lorentzian thus allowing 
I P'-d p I to stay bounded for finite P'. Such an imaginary parts does not arise as long as the collision 
term is neglected. 

The physical meaning of the variable P can be clarified by analyzing the transition amplitude between 
two states, < | and | g^,2 >, that contain at least one, reducible, on-shell particle with initial 

incoming momentum p, and final outgoing momentum q. In our convention the momenta are directed 
towards the interaction vertex. By definition 

p" - {-q"} = pi". (2.55) 

So P is the momentum loss due to interactions with a non-trivial background. The fact that for homo- 
geneous systems P = indicates that the interpretation is of a statistical nature. In homogeneous systems 
the particles will interact with the background, but the effects will average out. It is easy to show that 

2 2 -2 2 . u u /o r:a\ 

p = q = p — m p ~ q . (2.56) 

If G'^^{P, p) is restricted to being on-shell, with p on-shell as well, then it contains no dynamical infor- 
mation on the inhomogeneities. 

Thus far, we have derived the kinetic field theory for self interacting bosons. In the homogeneous 
limit the functions G°'^{P, p) describe the distribution of on-shell bosons over the available energy and 
momentum states. In this limit, p is the momentum of the propagating bosons. When we only require 
the asymptotic states to be on-shell P can be interpreted as the loss of momentum due to interactions 
with inhomogeneities in the system. 

In the next section, we will discuss a symmetry property of non-equilibrium field theories. It generates 
identities among Green functions that, among other things, express the impact of initial conditions on 
these Green functions. It serves to yield a satisfying interpretation of the 2-point functions and their 
arguments beyond the limit of quasi-homogeneity. 



12 



3 The Broken SO (1,1) Symmetry of Non-Equilibrium Field The- 
ory 



3.1 SO(l,l) Symmetry 

In this section, we focus on a symmetry underlying the non-equihbrium dynamics of the scalar field. But 
first let us make some general remarks. The transition from standard quantum field theory to the CTP 
formalism gives rise to a doubling of the number of fields. Every physical field becomes a CTP doublet, 
and on the internal space of these CTP doublets, we have defined the degenerate metric Cab- Now consider 
a theory containing an n-tuplc of fermions (and/or bosons) symmetric under SU(n) transformations. The 
free CTP action will then, through doubling, be symmetric under SU(n,n). As expected, this symmetry is 
broken, although an SU{n) SU{n) subgroup remains. Now SU(n,n) has (2n)^ — 1 generators, whereas 
SU{n)^SU{n) only has 2?!^ — 2 generators. Consequently the Noether theorem will yield 2rt^ + 1 
currents whose non-conservation can be computed exactly. For a broken symmetry, Ward identities 
can be derived and considered useful if, either the symmetry breaking terms are small or, symmetry 
violation is calculable. Here we will consider the symmetry relevant to the Xip'^ model. The real scalar 
field has a trivial S0(1) symmetry, which under CTP doubling becomes a broken S0(1,1) symmetry. To 
our knowledge this has gone unnoticed, although there is a relation to the SU(1,1) lie algebra used in 



quantum optics |29 . This section focusses on the S0(1,1) symmetry of the free real scalar field. We will 
give the Noether current, the Ward-Takahashi identities, discuss its breaking by interactions and initial 
conditions. 

The CTP action for the free field 

^CTP[^1 - I d*x ^Cabid^W^'d^^ip'' - mW'} (3.1) 

is invariant under global S0(1,1) transformations. The generator of these is 

and an infinitesimal transformation applied to the fields (p" yields the change 

Sip" = arfcV'' , I a |« 1 , (3.3) 
The associated S0(1,1) Noether current is 

r = eab^°9V' ■ (3.4) 

with tab the completely antisymmetric two by two CTP tensor. Note that (Pxj^{x) is a dimensionless 
quantity. In the general case, the expectation value of the Fourier transformed current is given by 

3''{P) - ^ P''{G^\P, P) - G^'iP, p)} . (3.5) 
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Its zero'th component is related to the canonical commutator by 



fix,t) = limj [ ^{x,t) , 7r{y,t) ] ) , (3.6) 



SO that we have 

j'iP) - J ^J'(P) (3-7) 
from canonical quantization. It implies that j'^{X) is an infinite constant. To see the meaning of this, 
consider the expectation value of the commutator of two fields 



mx),ip{y)])=G'\x,y)-G''ix,y). (3.8) 

It represents a pulse arriving at the point x — y from the backward lightcone, and subsequently going 
out into the forward lightcone ]30| . The temporal charge density is the amplitude corresponding to this, 
summed over all possible momenta. Put differently, j'^{x) is the amplitude density for absorbing a particle 
at X and subsequently emiting it. That it is infinite is due to the fact that it is obtained from 

-±^{G^\X, x)-G'\X, x)} , (3.9) 
a x^ 

evaluated at xf' — 0. It represents the amplitude density for absorbing a particle and reemmiting it at a 
relative distance x. Physically, at xf = the only contribution can come from 5 = 0, yielding the delta 
function. The infinity is a result of the expansion of these amplitudes in terms of plane waves. The use 
of extended wavepackets on the other hand would yield finite expressions. We see that the propagation 
of particles is represented by subsequent absorptions and reemmisions. In the free field vacuum this 
amplitude is conserved. The temporal charge density is a creation/ annihilation density, whose free field 
vacuum value describes the free propagation of particles. It is sensible to make a vacuum subtraction 
setting the free field temporal charge equal to zero, which corresponds to normal ordering. 

Analyzing the eigenstates of the S0(1,1) generator sheds light on the interpretation of j'~'{P) in terms 
of particles. The fields 

are normalised eigenvectors of the S0(1,1) generator. These are the generators of the Keldysh base i.e. 
the physical representation. [|l^. We assign the eigenvalues q± — ±1, or so-called temporal charge, to 
these. With respect to particle/anti-particle exchange, the states generated by ip~^ are symmetric and 
those generated by ip~ are anti-symmetric. Since, for neutral bosons, particles and anti-particles are 
indistiguishable, the 1-particle, negative temporal charge states are redundant. In the physical limit, the 
field carrying negative temporal charge < if>~ > vanishes. The field-operator ip'^ creates a symmetric 
state out of the vacuum, representing a physical particle. 

Let us see if the propagator is consistent with the demand that physical particles allways move into 
the forward light-cone. The propagator is given by 7p{(p^ (p^}. By considering the retarded and advanced 
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2-point Green functions, given by 

Gret{x,y) = {ip+ {x)ip- (y)) ^^^^^ 
Gadv{x,y) = {lp+ {y)ip^ [x)) 

we can identify the propagators for the particles with + and - temporal charge respectively. They 
represent respectively the outgoing and the incoming particles. In terms of these functions we can write 

f{P) = / ^ r{G.e*(P, P) - Ga,4P, P)} (3.12) 

showing that the temporal current is directly related to the absorptive part of the 2-point Green function 
Thus we find that 

{G^HP; p)-G'^{P, p)} (3.13) 

is the creation/annihilation density, in (P, p)-space, corresponding to momentum-loss P and the momentum- 
flow p. Since we have lost all coordinate space information, we cannot say in what region of coordinate 
space the momentum flows and where it is lost. However such questions can be addressed by constructing 



appropriate wave-packets. The commutation relation Eq.(3.6) expresses that if one integrates out the 
energy-flow, then P = 0. Put differently, the momentum is lost to other particles with different energies 
but it is not destroyed. The temporal current is the total flow of the creation/ annihilation density in 
terms of either a location, j'^{X), or the energy-momentum dissipation j'^{P). In detailed balance, the 
vacuum subtracted flow will be zero. The total temporal charge 



Q{t) = j d^x eab(^"5V' , (3.14) 

at time t is the total creation/ annihilation amplitude. Its Fourier transform is 

Q(e) = J dte'^'Qit) (3.15) 

the amplitude associated with the energy-loss e. The divergence of the temporal current is thus the 
amplitude for the energy-momentum loss P due to all interactions. 

The divergence of the temporal current is related to the integrated kinetic equations for G^^ and G^^ 
since 

P-^fiP) = / ^ P.P-'iG'HP, P) - G'HP, p)} . (3.16) 
As expected, and as can be seen from Eq.( ^.51 ), the global S0(1,1) symmetry is broken by the Xip'^ self- 



interaction. There is, however, another source of S0(1,1) symmetry breaking. Consider the free-field, 
vacuum, 2-point functions. If S0(1,1) symmetry were realized then G"^'' would be expressible as 

<,„(p)=5ic'^' + 52e'^'' ■ (3.17) 
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with gi and g2 functions of the momentum p. The actual Green function can be decomposed however 
to display the form 

= ( p2 _ ^2 _ ^Q^-i^ab _ 2m9{ f)6{ p2 - TO2)e"^ - 27n5{ p'' - m'') I ^ ^ ] (3.18) 

where one sees that the third term explicitly breaks SO (1,1) symmetry. A proper definition of the 
non-equilibrium theory, i.e. fixing the initial conditions, also requires the symmetry to be broken. 

3.2 Initial Conditions Breaking SO(l,l) 

For the sake of simplicity, we examine the symmetry breaking by a 2-point kernel Kab{x,y). This covers 
both equilibrium and non-equilibrium initial conditions, that describe an inhomogeneous system of free 
bosons. If A = 0, then the divergence of the temporal current is 

d^f{x) = J d^y K,b{x,y)^\y)Tl^\x) . (3.19) 

If this is integrated over x, after vacuum subtraction, we obtain the the total particle production ampli- 
tude, 

A^tot = i^a6o{a"^-a:L}. (3.20) 

These particles are created to form the initial state. Taking the most general form for the symmetry- 
breaking kernel to be 

Kab{x, y) = ki{x, y)dab + k2fab (3.21) 



with 



and using 



fab = r^, (3.22) 



we see the divergence of the temporal current is directly related to the 2-point correlation function 

duf{x) = j d^y{ki{x,y) + k2{x,y)}Gcorr{x,y)} . (3.24) 

In terms of temporal charge fields we have 

Gcorr{x,y) = {ip+{x)ip+{y)) , (3.25) 

representing the correlation among the particle production at x and at y. So far wc have discussed the 
expectation values of the S0(1,1) current operator. However in the context of Quantum Field Theory 
symmetries, even explicitly broken symmetries, generate a whole hierarchy of identities among Green 
functions. 
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By applying the usual techniques to the functional integral, assuming the absence of anomalous brea- 
king of the S0(1,1) symmetry, we can derive Ward-Takahashi identities for the broken symmetry. In the 
case of initial conditions symmetry-breaking (ICSB) discussed above we find 

^ n n n 

/n 
d'y Kab {z, y) {ip'^ (z)^" (y) J] ^'^^ {x,)) , 

(3.26) 

where the braces indicate functional integral averages. For n = we recover Eq.( 3.1i^ . For n > we 



expect the first term to become small as z — > oo and the aapplication of Gauss theorem to the l.h.s. of 



Eq.(3.26) leads to the relation 



J2t^H^Hx,11^^^{x,)) = -X,,o . (3.27) 

For n=l,2 we have explicitly 

^'^ {x) = -Ktciy, z) o ^^™(y, z, x) , (3.28) 

and 

G,or{x,y)^~Kic{z,z')og''-^\z,z',x,y), (3.29) 

representing the contribution of the initial conditions to particle distributions at later times. 

We have seen how initial conditions break S0(1,1) symmetry in order to prepare the initial state, and 
we have calculated the effect this has on the CTP Green functions. The self-interaction also breaks the 
symmetry and we will deal with that next. 



3.3 Interactions Breaking SO(l,l) 

The divergence in the temporal current due to interactions is 



d.f{x) = ~j^habcdTtv''{x)^\x)^^{x)^^{x) . (3.30) 

This explicit symmetry breaking (ESB) is local and therefore persistent, unlike the case discussed above. 
A simple example of the violation of particle number by interactions is the decay of one boson into three 



bosons. Note that Eq.(3.30) is just the d p integrated kinetic equation for G — G . The integrated 



Vlasov- and collision terms can all be recovered from Eq.( 3.30 ), and interpreted along the lines of (off- 
shell) particle creation and annihilation in the vacuum of the interacting theory. 
For ESB we also give a Ward-Takahashi identity, assuming no ICSB one finds 

n n n 

^ (r(^)n^"^(^:')) = Y.^^\z~x,)T:^^\z)^p^^^(^xi)) 

A " 

- ^WT,^(^"(^)^''w^^w^^wn^"^(^j-)) ■ (3.31) 
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Once again for n > 0, we expect an integrated identity to hold. We find 



A 

^^a,ga,...a,_,ca, + ^...a^ ^ ^/iabcd O T.^'^^'^''"""! ■ ■ . (3.32) 
J = l 



Again we give to explicit examples, for n = 1 



^'^ ^^hbcdeOTfT'gg''''3f , (3.33) 

and for n = 2, we find an equation for the correlation function 

Gcorrix, y) = J^ habcdrt o ^'^"-^ . (3.34) 

In any realistic case, both mechanisms of S0(1,1) breaking contribute. The broken S0(1,1) symmetry 
generates a hierarchy of equations for the Green functions. The generalization of Eq.( |3.26 ) and Eq.(3.31) 
to higher-order kernels Kabc... and other polynomial interactions is straightforward. 
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4 Energy-Momentum Conservation 

In this section, we show that the divergence of the canonical energy-momentum tensor vanishes due 
to the kinetic and constraint equations. It wiU turn out that the true difference between Vlasov- and 
coUision-terms Ues in the (non-)locahty of the interactions that they represent. 
The energy-momentum tensor in the Xip^ model is given by 

= To'^ + ^TZ, (4.1) 

where Tq'^ refers to the free field energy- momentum tensor and to the interacting part. In terms of 
2-point fimctions this can be written as 

(^o^'') - ^ / -^i-lP'P" + - l^'^i f - - ^^)}Gcorr{P, P) , (4.2) 

where Gcorr has been defined in Eq.(3.23). In deriving this expression, linear terms in p'^ have been 
dropped, because they do not contribute to the integral since Gcorr{P, p) is symmetric in p'^. The 
interaction contribution to the energy momentum tensor is connected with the coincidence limit of the 
four point Green functions. One possible choice is 

(O = lv''-g''^Hx,x,x,x)+g^'^\x,x,x,x) . (4.3) 

It is important to recognize that the interaction part also contains the contribution of the independent 
variable oc^^,^^- Energy-momentum conservation requires 

p,{Tn^^,p,{Tt:t) ■ (4.4) 

We will investigate this relation for three cases. First we will consider the free field situation, followed by 



a verification of Eq.(4.4) taking into account the Vlasov interaction and, finally, treating the full theory. 



Our procedure will be to use the kinetic and constraint equations to calculate the l.h.s. of Eq.(4.4) from 



Eq.(12) and conclude the equality with the r.h.s. of Eq.(4.4) 



For the free field case, A = 0, the kinetic and constraint equations read 

p.PC'^P, p) = 



(4.5) 

{ p2 + lp2 _ ^2|gafc(p^ ^ c''\2nfS^{P) . 



Applying to Eq.( |4.2| ), it is straightforward to see that the l.h.s. of Eq.(4.4) vanishes. In the Vlasov 
approximation, we neglect the collision terms, i.e. we treat the 4-point Green function as a disconnected 
product of 2-point fimctions. So the interaction term is 

(TZ) = ^v''''{G'Hx,x)G'Hx,x)+G^\x,x)G'\x,x)} . (4.6) 
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The factor 3 counts the possible ways of writing a 4-point Green function as a disconnected product of 
2-point Green functions. The kinetic and constraint equations are now given as 

^{ + kP' - m^}G^'{P, p) - I / ^'^G--[P\ p') cosh{\P' .d p}G'^\P - P\ p) (4.7) 
= ic''\2tt)^5^{P) 

Applying to Eq.( [4.2| ) and using the Vlasov equations yields 

'4 p' ^4 ^1 



1 /■ p f (PP' (p p 



(27r)4 J (27r)4 (27r)4 



pn G^'{P', p')smh{-P'.dp}G'HP~P\ P) 
+ G'\P', p')siTih{\p'.dp}G^\P~P', p)} 



+ Pl{ G^HP', p')cosh{^P'.dp}G'\P^P', p) 



G^\P', p')cosh{^P'.dp}G^^{P-P', p)} . 



(4.8) 



The p-integration in the terms containing the sinh and cosh functions can be performed. For the sinh 
terms this yields 



"^'^ p'^G'^'^iP', p') smh{l-P'.dp}G'^\P-P', p) = -i^G-(P', p') 



(2^) 



P ^a, 



{2n) 



G'^\P-P', p), 
(4.9) 



to be integrated over P' and p'. For the cosh terms one can use 

^^G^^iP', p') cosh{\p'.dp}G'^\P-P\ p) = -:^G-(P', p') 



P na. 



G'^'iP-P', p). 

(4.10) 



These terms generate the r.h.s. of Eq.(4.4). Note that, since the Vlasov term has gradient operators 
between a p-independent factor and G''^(P, p), the gradients play no role in energy-momentum conser- 
vation. 



In the fully interacting case, the confirmation of Eq.(4.4) becomes difficult as long as we stay in 



(P, p)-space. Substitution of Eq.(£^5l|) and Eq.( |2.52| ) into the l.h.s. of Eq. now yields 

P^(rr(P)) ^HJ$^.i^.P''{^\{P-P'. P)A-G'^HP', p) + E^(P-P', p)A-G^HP', p) 
-G\iP~P', p)A+S=i(P', p)-G\iP~P', p)A+S=2(p/^ 
^{I]2^(P-P', p)A-G'\P', p) + Si,(P-P', p)A-G'^iP\ p) 
+G\{P - P', p)A+S=i(P', p) + G\{P - P', p)A+S=2(p/^ (4 

where we have allready used the fact that the free field and Vlasov terms either vanish or are taken care 



of by corresponding terms in the interaction part. The d p-integration in Eq.(4.11) is not trivial, as it 



20 



was for the Vlasov term. All the gradients must be included. We do this by transforming Eq.(4.11) back 
to coordinate-space. There we have 

+ {^[G^,{y,z)^^\z,x) + G\{y,z)^'^\z,x)]} , (4.12) 

which we can identify as the point-split implementation of the product rule by using the symmetry 
of the Green functions under the coherent exchange of CTP indices and the corresponding space-time 
arguments. Consequently we have the divergence of the free field energy- momentum tensor expressed in 
terms of the connected 4-point Green function 

d,{T>;-) ^ ^d'^{G^''H^,^,^,^) +G'''\x,x,x,x)} . (4.13) 



Using Eq.(4.6) we see that Eq.(4.4) is satisfied, so that in general energy and momentum are conserved. 
Contrary to energy- momentum conservation for the Vlasov terms we required all the gradients. The 
physical reason for this lies in the nature of the interactions that these terms represent. The Vlasov-term 
results from local interactions with (off-shell) particles from the background. The interaction vertex 
is the point-like Xip"^ interaction. The collision terms contain non-local interactions, like, for example, 
the two-boson exchange interaction. The gradients probe the locality of these interactions. Neglecting 
gradients thus violates energy-momentum conservation, unless specific non-local interactions are added 
to the Lagrangian density. 

Energy-momentum conservation is guaranteed by the kinetic- and constraint equations. This points 
to the fact that any inhomogeneities occurring in such a field theoretic model come from a spontaneous 
breaking of translation invariance. A gradient-expansion violates locality if used in the collision term. A 
Vlasov-approximation is restricted to the local response of the system to inhomogeneities. The temporal 
current non-conservation is a differential equation for the difference G^^ — G^^ in the coincidence limit and 
energy-momentum conservation is a differential equation for the sum G^^ + G^^ in the coincidence limit. 
The two equations, given apropriate boundary conditions, consequently determine the coincidence limit of 
the 2-point Green functions. Put differently, energy-momentum conservation and temporal current non- 
conservation determine the U.V. behaviour of the dressed propagators. In the next section we analyse a 
free field case in the presence of inhomogeneities. 



5 Free field 

In the previous sections we have seen that energy-momentum conservation and temporal current non- 
conservation equations determine the coincidence limits of G^^(a;, y) and G^^(x, y). Inhomogeneities have 
played a major role in establishing the conservation since they reflect the energy-momentum exchange 
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between the particles. Inhomogeneities on the other hand make their appearance in the gradient expan- 
sions of the interaction terms in the kinetic and constraint equations, but also they appear quadratically 
in the free field equations. In this section we will concentrate on the free field equations 



(5.1) 

To find the solution to Eq.( p.l[ ) for arbitrary initial data, we exploit the linearity of these equations. It is 
used to discuss the effect of inhomogeneities and off-shellness. Following this we give an explicit example 
of an inhomogeneous non-equilibrium system. 

5.1 General solution 



By inspection of Eq.( |5.1[ ) one sees that the decomposition of G"'' into a vacuum contribution, G^J^^, and 
a medium term, G"^, 

G^\P, p) = GtAP. P) + G^\P. P) , (5.2) 
solves the free field equations, where G°|^^ is 

GtAP. P) = (27r)^c°^ _^ //^P ^ (5.3) 

+ jP^ — TO^ 

and G'^''{P, p) satisfies the set of homogeneous equations 

p.PG-\P, p)^0 

(5-4) 

{ p2 + lp2 _ TO2|(^afc(p^ . 

Since G°|j^ describes the propagation of free particles, their P-dependence is a simple 5^{P) as ex- 
pected from the absence of interactions. All the momentum- loss information is included in G"''. It is 
the momentum-loss suffered by the particles during the preparation of the initial condition. Particle 
propagation at later times will not couple to these initial inhomogeneities and thus not dissipate them. 

Let us consider that G"''(X, x) is given by an initial value at X° = and x'^ — and denote it as 
Gq''{X^, x^). Since we work in momentum space, it is useful to introduce the Fourier transform of the 
initial conditions Gg^(P*, p^) 

Gf{P\ f) ^ J^X J xe^'P^^'-' Po ^'Gg''(X^ x') . (5.5) 



This Fourier transform can also be decomposed in the fashion of Eq.(||2|) as G^^ = -f Gf'. The 

linearity of the equations implies the solution for G"'' has the form 

&\P, p) = Gl\P, p)5{ p.PnSi P' + \P^ - m") . (5.6) 

Due to the symmetry, G^^{P, p) — G^^{P, — p), of the off-diagonal 2-point functions and the relation 

G"(X, x) = 0{ x°)G21(X, x) + 9{- i°)Gi2(X, x) , (5.7) 
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specifying the diagonal elements in terms of the off-diagonal elements, the four functions G contain 



only one unknown function g{P^ p), satisfying Eq.(5.4). In thermal equilibrium, it is well-knowm that 
g{P, p) is given as 

5(P, p) = - m')6\P) , (5.8) 

where gBE{p^) is the Bose-Einstein distribution function. 

Consider the Wigner representation of the general solution 

5(X, p)= j i!^e^^-^^<S( p^P>^)5{ p' + _ ^2)^(p^ -) (5 9) 

Since we deal with massive particles p^ ^ and two of the four P integrations, say over P" and the 
component parallel to p> 

P|f = ^,f , (5.10) 



P^ P 



can be performed using 



5{p^ + \P^~m^) = ^{<5(Pt|-27a) + ^(P||+27a)} 
where we have defined the measure of off-shellness a by 



(5.11) 



and 



p2-,n2-ip2 (5.12) 



(5.13) 



1 Pi P" 



Pi 



and used the first identity in Eq.( 5.11 ) to obtain the second. In these coordinates p has only two non- 
vanishing elements; p^ and p. Thus 7 is a Lorentz boost factor and a is not a Lorentz- invariant quantity, 
but rather a frequency. After performing the integrals in P'^ and P|| , the integrand will have the general 
form 

9{X, p, p") ^ J j^—^W{p, p'',-f,a,X±_)e P +W{ p, p'\-f,a,Xj_)e 7^ 

(5.14) 

Eq.(5.14) is a superposition of waves with frequencies 7a moving along the X||-axis. The Lorentz boost 
factor relates these waves to spatial oscillations in the rest-frame defined by p = 0. 

Although there is no interaction among the particles, there are collective off-shell effects due to the 
preparation of the initial condition. In a free field theory, particles with a pure momentum are described by 



plane waves. Eq.(5.14) tells us that a particle cannot be in a pure momentum state if there are transversal 
inhomogencities (T.I.'s). Clearly they would disturb any plane wave into a more complex wave. From 
the reality condition on a, we see that T.I.'s must vanish for p to be on-shell. The constraint p^P^ = 
implies that the rest-mass of the particles is unaffected by inhomogencities. The energy-loss is solely due 
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to longitudinal inhomogeneities (L.I.'s), in contrast to the T.I.'s, responsible for removing the particles off 
the mass shell. By going to the rest frame p = of the particles moving along X|| , one sees that > m 
for non- vanishing T.I.'s. This implies that the particles will not reach asymptotic on-shell states unless 
they were not deflected in the first place. One of the two remaining integrations over the components of 
Pi can be replaced by an integration over the off-shellness. Integrating out the a represents a summing 
over quantum fluctuations dressing the particles. 



5.2 Free expansion of a Bose-Einstein gas 

In this subsection, we give an example of a system as discussed above. The free expansion of a Bose- 
Einstein gas will be treated. Consider (3 + 1) dimensional Minkowski space divided into two half-spaces, 
and Af_, by a Dirichlet-wall at ~ 0. The gas in Af_ is assumed to be in thermal equilibrium at a 
temperature T, while in Af+ a vacuum persists. The total-space propagator G{x,y) is defined as a sum 
of two half-space propagators 

G{x,y)^e{x^)e{y^)G+{x,y)+e{-x^)e{-y'')G-{x,y) (5.15) 

where G+ is a vacuum propagator and G_ is a thermal propagator. These can be computed from total- 
space vacuum and thermal propagators by the method of images. This procedure can be justified by an 
eigenfunction expansion in terms of solutions of the half-space Klein-Gordon equation. For the half-space 
vacuum, one can express the associated Green function in terms of the total-space vacuum propagator 
Gvac{x, y) 

G+{x, y) ^ G^ac{x; y) - G^ac{x\ -x^; y) . (5.16) 

By using the fact that the reflection {x^,y^) {—x^,y'^), is equivalent to {X^, x^) — > (— ^ x^ ,—2X^), 
in the centre of mass (cm.) and relative coordinates, we can write down the relation in the Fourier 
representation as 

G+{P; p)^G.,ac{P; p)-G,„ac(P\-2p3; p\-^P^) . (5.17) 
In the thermal half-space on the other hand we have, using the thermal total-space propagator Gt^P, p), 

G.{P- p)^Gt{P; p)-Gt{P\-2p'; ■ (5.18) 

. These propagators describe a stationary state which at — x'^ — coincides with our initial condition. 
In momentum space the initial conditions in each of the half-spaces is now represented by the apropriate 
propagator integrated over P° and Thus we compute the initial condition for each half-space 

r dP° d f)° 

G± o{P\ P') = / ^^G±(P, p) . (5.19) 



The total space initial condition is calculated from the functions G± o according to Eq.(5.15). The 
product with the step functions in coordinate space becomes a convolution in momentum space. The 
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Fourier transform of the product of step functions is found by noting that 

0{±x^)9{±y^) = e{±X^){d{±2X^ - x^) - 9{t2X^ - x^)} , (5.20) 

and Fourier transforming this yields 

_ _ ±4 p3 

^ [ p3 + ie][(p3)2 „ 4( p3)2 ± j2P3e _ £2] ' (5-21) 

where the hmit e, e — > is understood. Thus the required fuU-space initial conditions are given by 
Go{P\ f) = I ^^h{P', P')G+ o{P^ - Pi P, - P,) 

J 111 2,71 



^^^e^iP', p')G^ o(P. - Pi P, - P,) ■ (5.22) 

111 2,71 J J 

Having now determined the initial condition Green function, we focus on the calculation of the medium. 
Using the linearity of the free field equations and the Fourier transform we can subtract all vacuum 
contributions from Go{P'', p') to obtain Go{P'', p'). By making a high-tcmperature expansion 

9BEiio) ~ , iPio) « 1, (5.23) 
pui 



these integrals can be calculated analytically. Repeating the steps leading from Eq.( |5.9| ) to Eq.(5.14) yields 
the general solution. In the case under consideration inhomogeneities only exist in the X'^ direction, so 
the dependence on P± is simply a S{Pj_). The measure of off-shellness now takes the form 

a = a/ p2 - , (5.24) 

and the boostfactor is 

(5.25) 



_ Pt 

'pI 

The solution thus has a simpler form 

g{X, p) = ^^^{(A)sin7a(X3-^X") + (i?)cos7a(X3-^X0)} , (5.26) 
70; p p p 

where the A and B are given by 



—4 9 9 —2 

7" Pq - 7 a Pi 



, pI- pI - il Oi 
B = -ttT^^^j ^^^2 • (5-27) 

Po - 7 a Pi 



If one averages Eq.(5.26) over small segments in X^, or in p^ , the a = poles become stronger by 
a factor a^^. This suggests that the semi-classical limit may be uncovered not by setting a = 0, but 
rather by some kind of averaging. Since the propagation of the waves appears to be dominated by the 
pole at a = 0, this suggests averaging over the off-shellness. In the general case this could replace an 
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integration over P±, but here there is no non-trivial P± dependence. Consequently we will integrate over 
the momenta perpendicular to the axis when averaging over a. The integration over a is restricted 
to the range 



< a < y P5 - pi-m^ (5.28) 

due to the reality condition on a. Using the freedom to make a Lorentz transform to the p'^ frame, one 
obtains 



^3, Pq) = 47r— ^ 

Po Po 
1 



2j_i(^^, p,X,) 

P" } , (5.29) 



[2 J Pl - m2X3COs(2 J pl - m^Xs) + sin(2 J _ m2x^)] 



X3PI 

where we have defined the function j_^(a;, y) by 



j-uix,y) = I dz {-yT{\/l- z^y-i cos(yz) . 



(5.30) 







with ly = 1. The function .9(^3, Pq) is related to the non-equilibrium generalization f{X^, p^) of the 
Bose-Einstein distribution gB e by 

/(X3, Po) = plgiX,, Po) , (5.31) 

except for an overall normalization. The distribution function /(X3, Pq) is represented graphically in 
fig. 2 as a function of X3 for fixed fP, and in fig. 3 as a function of p'^ for fixed X3. Two interesting 
limiting cases that one may consider are po ^ 00, and X3 ~* 0. At increasing energy, /(X3, Pq), 
behaves more and more like a step function being zero in M_|_ and thermal in Af_. This is the result 
that one expects classically. The second limit, X3 0, represents a purely field theoretical effect. The 
divergence near the Dirichlet wall is nothing but the one-plate Casimir effect. In the presence of such a 
wall, quantum fluctuation near it will be modified. The modification of vacuum fluctuations gives rise to 
a 1/^3 behaviour [^ . On dimensional grounds one expects the finite temperature contribution to be 
proportional to T/X3, which is the behaviour we find. An explicit calculation verifies this interpretation. 
As the wall is instantaneously removed at t = 0, this leads to particle creation. It is exactly this type of 
effect that the S0(1,1) Ward identities for ICSB adress. 



6 Conclusions 

In principle a kinetic field theory of non-equilibrium physics is viable beyond a two-scale approach. It is the 
ideal starting point for treating inhomogeneous systems out of equilibrium. In particular we have studied 
a transport theory for the scalar Xip^ field theory, constructed using the closed time path formalism. 
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In this case we find that energy momentum conservation is guaranteed by the dynamics of the 2-point 
functions. Gradient expansions occurring naturally within the formalism cannot be truncated without 
either violating energy-momentum conservation, or introducing non-locality. It has been shown how the 
broken S0(1,1) symmetry of non-equilibrium field theory is useful in establishing a physical interpretation 
for G""^. Temporal current non-conservation combined with energy- momentum conservation determines 
the ultra-violet behaviour of the dressed propagators. The free field equations have been analyzed on 
their non-equilibrium content. Inhomogeneity and off-shellness are intimately related in this context. 
This is not expected to change in the interacting theory. By integrating out the off-shell contributions 
one arrives at quasi classical distribution functions. 

Acknowledgements 

Wc arc grateful to J. Hiifncr for his continuing encouragement and discussions during the course of 
this work. One of us, F.W, would like to thank C. Villarreal for useful discussions on the Casimir effect 
at finite temperature. 



27 



References 

[1] L. P. Csernai, J. I. Kapusta, Gy. Kluge and E. E. Zabrodin, Proc. of Quark Matter '91, Nucl.Phys.A 
544, (1992). 

[2] Blaettel et. al, Rep. Prog. Phys 50, 63 (1993). 

[3] A. G. Cohen, D. B. Kaplan and A. E. Nelson, Phys. Lett. B 336 , 41 (1994); J. M. Cline, Univ. of 
Minnesoty preprint, UMN-TH-1259-94, (1994). 

[4] J. D. Bjorken, Stanford preprint, SLAC-PUB-6488, (1994). 

[5] E. Calzetta and B. L. Hu, Phys. Rev. D 37, 2878 (1988). 

[6] L. V. Keldysh, JETP 20, 1018 (1965). 

[7] J. Schwinger, J. Math. Phys. 2, 407 (1961). 

[8] R. Feynman and F. L. Vernon, Ann. Phys. 24, 118 (1963). 

[9] H. Umezawa, H. Matsumoto, and M. Tachiki, Thermo Field Dynamics and Condensed States (North 
Holland, Amsterdam, 1982). 

[10] A. J. Niemi and G. W. Semenoff, Ann. Phys. 152, 105 (1984). 

[11] L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974). 

[12] K. Chou, Z. Su, B. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). 

[13] J. M. Eisenberg and G. Kaelbermann, Phys. Rev. D 37, 1197 (1988). 

[14] P. Carruthers and F. Zachariasen, Phys. Rev. D 13, 950 (1976). 

[15] B. L. Hu, Physica A 158, 399 (1989). 

[16] H. Klcincrt, Forts. Phys, 30, 187 (1982). 

[17] P. Danielewicz, Ann. Phys. 197, 154 (1990); P. Daniclcwicz, Ann. Phys. 152, 239 (1984). 
[18] H. Th. Elze and U. Heinz, Phys. Rep. 183, 81 (1989) ; U. Heinz, preprint NSF-ITP-94-03. 
[19] Y. Zhang and L. Willets, Phys. Rev. C 45, 1900 (1992). 
[20] A. Abada and J. Aichehn, University of Nantes preprint, LPN 94-02. 
[21] W. Florkowski, University of Heidelberg preprint, HD-TVP-94-12. 
[22; 



F. Cooper, Y. Kluger, E. Motola and J. P. Paz, |hep-ph/9404357 



28 



[23] W. Botermans and R. Malfliet, Phys. Rep. 198, 115 (1990). 

[24] R. Feynman, Statistical Mechanics; A Set of Lectures, (W. A. Benjamin Inc., 1972). 
[25] A. Dannenberg, Phys. Lett. B 202, 110 (1988). 

[26] F. M. C. Witte and M. G. Schmidt, University of Heidelberg preprint, HD-THEP-94-44. 

[27] B. L. Hu and H. E. Kandrup, Phys. Rev. D 35, 1776 (1987) ; H. E. Kandrup, J. Math. Phys. 28, 
1398 (1987); H. E. Kandrup, Physica A 158, 336 (1989) 

[28] H. A. Weldon, Phys. Rev. D 28, 2007 (1983) 

[29] M. Ban, Phys. Rev. A 47, 5093 (1993); M. Ban, Phys. Lett 193, 121 (1994) 

[30] P. A. M. Dirac, Lectures on quantum field theory, Belfer Graduate School of Science, Yeshiva Uni- 
versity, N.Y. (1966). 

[31] S. A. Fuller, Aspects of Quantum Field Theory in Curved Spacetime, Cambridge University Press, 
(1989). 

[32] M. Veltman, Diagrammatica, Cambridge University Press (1994). 



29 



Figure Captions 



• Fig. la and lb: The diagrammatical representations of Eq.(2.24) and Eq.(2.25) are given. The 
shaded triangle denotes the IPI 3- vertex qJ^^^, an open triangle refers to the general 3-point Green 
function C?"*"', thin solid lines represent outgoing lines, single blobs correspond to the mean fields 

c^" , and lines with a blob along the line represent the full propagator G""^. Additionally in lb the 
shaded disc is the IPI 4- vertex a^;,^^ and the open disc the general Green function Q"-^'^'^. 

• Figures 2 a-d: The distribution function f{x, e) as a function of the energy e = pQ^tx = —5, —1,1,5 
fm/c. Here the temperature T = 2 GeV and the mass of the bosons was set m = 0.2 GeV. 

• Figures 3 a-d: The function f{x, e) as a function of a; at e = 0.2, 0.3, 1, 2 GeV. 
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